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Abstract 

Measurements of 6-hadron masses are performed with the exclusive decay modes 
B+ J/^K+, J/ipK*^ , J/i^K^, 5° ^ J/ip(j) and ^ J/V'A 

using an integrated luminosity of 35 pb~^ collected in pp collisions at a centre-of- 
mass energy of 7TeV by the LHCb experiment. The momentum scale is calibrated 
with J/ip ^ IJ'~^IJ'~ decays and verified to be known to a relative precision of 2 
xlO~^ using other two-body decays. The results are more precise than previous 
measurements, particularly in the case of the Bg and A|] masses. 
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1 Introduction 



Within the Standard Model of particle physics, mesons and baryons are colourless objects 
composed of quarks and gluons. These systems are bound through the strong interac- 
tion, described by quantum chromodynamics (QCD). A basic property of hadrons that 
can be compared to theoretical predictions is their masses. The most recent theoretical 
predictions based on lattice QCD calculations can be found in Refs. [1, 2]. The current 
experimental knowledge of the 6-hadron masses as summarized in Ref. [3] is dominated 
by results from the CDF collaboration [4]. In this Letter precision measurements of the 
masses of the B~^, B^, and are presented as well as the mass splittings with re- 
spect to the B^. The results are based on a data sample of proton-proton collisions at 
y/s = 7TeV at the Large Hadron Collider collected by the LHCb experiment, correspond- 
ing to an integrated luminosity of 35 pb~^. 

The LHCb detector [5] is a forward spectrometer providing charged particle recon- 
struction in the pseudorapidity range 2 < t] < 5. The most important elements for the 
analysis presented here are precision tracking and excellent particle identification. The 
tracking system consists of a silicon strip vertex detector (VELO) surrounding the pp 
interaction region, a large area silicon strip detector located upstream of a dipole magnet 
with a bending power of about 4 Tm, and a combination of silicon strip detectors and 
straw drift-tubes placed downstream. The combined tracking system has a momentum 
resolution 6p/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c. Pion, kaon and 
proton separation is provided by two ring imaging Cherenkov (RICH) detectors whilst 
muons are identified by a muon system consisting of alternating layers of iron and multi- 
wire proportional chambers. 

The data used for this analysis were collected in 2010. The trigger system consists 
of two levels. The first stage is implemented in hardware and uses information from the 
calorimeters and the muon system. The second stage is implemented in software and runs 
on an event filter farm. Dedicated trigger lines collect events containing J/tp mesons. For 
this analysis all events are used regardless of which trigger line fired. 

Simulation samples used are based on the Pythia 6.4 generator [6] configured with 
the parameters detailed in Ref. [7]. QED final state radiative corrections are included 
using the Photos package [8]. The EvtGen [9] and Geant4 [10] packages are used to 
generate hadron decays and simulate interactions in the detector, respectively. 

The alignment of the tracking system, as well as the calibration of the momentum 
scale based on the J/ip ^ Z^'*'/^" mass peak, were carried out in seven time periods 
corresponding to different running conditions. The procedure takes into account the effects 
of QED radiative corrections which are important in the J/ip ^ /^^/^~ decay. Figure 1 
shows that the reconstructed J/ip mass after alignment and calibration is stable in time 
to better than 0.02% throughout the data-taking period. The validity of the momentum 
calibration has been checked using samples of — )■ tt'^7c~ , — )■ K~7t~^, — )■ K^n~ , 
4'{2S) — )■ /i"*"//", T(IS') — )■ /i"*"//" and T(2S') — )• decays. In each case the mass 

distribution is modelled taking into account the effect of radiative corrections, resolution 
and background, and the mean mass value extracted. To allow comparison between the 
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Figure 1: Reconstructed J/ip — )• /x"*"/!" fitted mass as a function of run number after the 
momentum calibration procedure discussed in the text. The vertical dashed lines indicate the 
boundaries of the seven calibration periods. A fit of a constant function (horizontal line) has a 
probability of 6%. The shaded area corresponds to the assigned uncertainty on the momentum 
scale of 0.02%. 
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Figure 2: Momentum scale bias q, extracted from the reconstructed mass of various two-body 
decays after the momentum calibration procedure described in the text. By construction one 
expects a = for the J/V' — )• M"*"/!" calibration mode. The black error bars represent the 
statistical uncertainty whilst the (yellow) shaded areas include contributions to the systematic 
error from the fitting procedure, the effect of QED radiative corrections and the uncertainty 
quoted by the PDG [3] on the mass of the decaying meson. The (red) dashed lines correspond 
to the assigned uncertainty on the momentum scale of 0.02%. 



decay modes, the deviation of the measured mass from the expected value [3] is converted 
into an estimate of the momentum scale bias, referred to as a. This is defined such that 
the measured mass is equal to the expected value if all particle momenta are multiplied 
by 1 — a. Figure 2 shows the resulting values of a. The deviation for the considered 
modes is ±0.02%, which is taken as the systematic uncertainty on the momentum scale. 
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2 Event selection 



A common strategy, aiming at high signal purity, is adopted for the reconstruction and 
selection of B+ J/^IjK+, 5° -> J/ipK*^, 5° ^ J/^K^, 5° ^ J/ipcj) and J/^A 
candidates (the inclusion of charge-conjugated modes is implied throughout). In general, 
only tracks traversing the whole spectrometer are used; however, since and A particles 
may decay outside of the VELO, pairs of tracks without VELO hits are also used to build 
Kg and A candidates. The P^r number of degrees of freedom (x^/ndf) of the track 
fit is required to be smaller than four. The Kullback-Leibler (KL) distance [ I 1] is used 
to identify pairs of reconstructed tracks that are very likely to arise from hits created by 
the same charged particle: if two reconstructed tracks have a symmetrized KL divergence 
less than 5000, only that with the higher fit quality is considered. 

J/ip ^ fi^fj^~ candidates are formed from pairs of oppositely-charged muons with a 
transverse momentum (px) larger than 0.5 GeV/c, originating from a common vertex with 
X^/ndf < 11, and satisfying |M^^ — Mj/^\ < 3a where M^^ is the reconstructed dimuon 
mass, Mj/^ is the known J/ip mass world average value [->], and a is the estimated 
event-by-event uncertainty on M^^. The selected J/tp candidates are then combined with 
one of K~^, K*^ — i- K^tt~, — K~^K~, Kg — t- tt+tt^ or A — > pir" to create 6-hadron 
candidates. Mass windows of iTOMeV/c^, ±12MeV/c2, ±12MeV/c2 (±21MeV/c2) and 
±6MeV/c^ (±6MeV/c^) around the world averages [->] are used to select the 0, Kg 
and A candidates formed from tracks with (without) VELO hits, respectively. Kaons are 
selected by cutting on the difference between the log-likelihoods of the kaon and pion 
hypotheses provided by the RICH detectors {AIiiCk-tt > 0). To eliminate background 
from — )■ J/ip(j) in the — t- J/ipK*^ channel, the pion from the K*^ candidate is 
required to be inconsistent with the kaon hypothesis (Aln£j<-_7r < 0). To further improve 
the signal purity, a requirement of pt > 1 GeV/ c is applied on the particle associated with 
the J /if) candidate. For final states including a {Kg or A), an additional requirement 
of L/ctl > 5 is made, where L is the distance between the 6-hadron and the decay 
vertex, and ctl is the uncertainty on this quantity. 

Each 6-hadron candidate is associated with the reconstructed pp primary interaction 
vertex with respect to which it has the smallest impact parameter significance, and this 
significance is required to be less than five. As there is a large combinatorial background 
due to particles originating directly from the pp primary vertex, only 6-hadron candidates 
with a reconstructed decay time greater than 0.3 ps are considered for subsequent analysis. 
A decay chain fit [12] is performed for each candidate, which constrains the reconstructed 
J/ip mass and, if applicable, the reconstructed Kg or A mass to their nominal values [■>]. 
The x^/iidf of the fit is required to be smaller than five. The mass of the 6-hadron 
candidate is obtained from this fit and its estimated uncertainty is required to be smaller 
than 20MeV/c2. 
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Table 1: Signal yields, mass values and mass resolutions obtained from the fits shown in Fig. 3 
together with the values corrected for the effect of QED radiative corrections as described in the 
text. The quoted uncertainties are statistical. 



Decay mode 


Yield 


Fitted mass 
[MeV/c2] 


Corrected mass 
[MeV/c2] 


Resolution 

[MeV/c2] 


B+ J/^K+ 
^ J/ijA 


11151 ± 115 
3308 ± 65 
1184 ± 38 
816 ± 30 
279 ± 19 


5279.24 ±0.11 
5279.47 ±0.17 
5279.58 ± 0.29 
5366.90 ± 0.28 
5619.19 ±0.70 


5279.38 ±0.11 
5279.58 ±0.17 
5279.58 ± 0.29 
5366.90 ± 0.28 
5619.19 ±0.70 


10.5 ±0.1 
7.7 ±0.2 
8.6 ±0.3 
7.0 ±0.3 
9.0 ±0.6 



3 Results 

The 6-hadron masses are determined by performing unbinned maximum likelihood fits 
to the invariant mass distributions, in which the signal and background components are 
described by a Gaussian and an exponential function, respectively. Alternative models 
for both the signal and background components are considered as part of the systematic 
studies. Figure 3 shows the invariant mass distributions and fits for the five modes 
considered in this study. The signal yields, mass values and resolutions resulting from the 
fits are given in Table 1. 

The presence of biases due to neglecting QED radiative corrections in the mass fits 
is studied using a simulation based on Photos [ ]. The fitted masses quoted in Table 1 
for the B'^ — )■ J /ipK^ and B^ J/ipK*^ are found to be underestimated by 0.14 ± 
0.01 MeV/c^ and 0.11 ± 0.01 MeV/c^, respectively, when radiative corrections are ignored; 
they are therefore corrected for these biases, and the uncertainty is propagated as a 
systematic effect. The bias for the B^ — )■ J/ip(f) mode is negligible due to the restricted 
phase space for the kaons from the (j) decay. There is no bias for the 5" — )■ J/ipK^ and 
A^ — 7- J/ipA modes since the J/ip, and A masses are constrained in the vertex fits. 

4 Systematic studies and checks 

To evaluate the systematic error, the complete analysis is repeated (including the track 
fit and the momentum scale calibration when needed), varying within their uncertainties 
the parameters to which the mass determination is sensitive. The observed changes in 
the central values of the fitted masses relative to the nominal results are then assigned as 
systematic uncertainties. 

The dominant source of uncertainty is the limited knowledge of the momentum scale. 
The mass fits are repeated with the momentum scale varied by ±0.02%. After the cal- 
ibration procedure a ±0.07% variation of the momentum scale remains as a function of 
the particle pseudorapidity rj. To first order the effect of this averages out in the mass 
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Figure 3: Invariant mass distributions for (a) J/il^K^ , (b) J/ijjK*^, (c) B^ — > 

J/tpKo^, (d) — )• J/i/j^, and (e) 5° — >• J/V"/" candidates. In each case tlie result of the fit 
described in the text is superimposed (solid line) together with the background component 
(dotted line). 



determination. The residual impact of this variation is evaluated by parameterizing the 
momentum scale as a function of i] and repeating the analysis. The amount of material 
traversed in the tracking system by a particle is known to 10% accuracy [L3]; the mag- 
nitude of the energy loss correction in the reconstruction is therefore varied by 10%. To 
ensure the detector alignment is well understood a further test is carried out: the horizon- 
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tal and vertical slopes of the tracks close to the interaction region, which are determined by 
measurements in the VELO, are changed by lxlO~^, corresponding to the precision with 
which the length scale along the beam axis is known [11]. Other uncertainties arise from 
the fit modelling: a double Gaussian function (with common mean) for the signal resolu- 
tion and/or a fiat background component are used instead of the nominal Gaussian and 
exponential functions. The effect of possible reflections due to particle mis-identification 
is small and can be neglected. Finally, a systematic uncertainty related to the evalua- 
tion of the effect of the radiative corrections is assigned. Tables 2 and 3 summarize the 
systematic uncertainties assigned on the measured masses and mass differences. 

The stability of the measured 6-hadron masses is studied by dividing the data samples 
according to the polarity of the spectrometer magnet, final state flavour (for modes where 
the final state is flavour specific), as well as whether the and A daughter particles 
have VELO hits. As a cross-check the analysis is repeated ignoring the hits from the 
tracking station before the magnet. This leads to an average shift in measured masses 
compatible with statistical fluctuations. In addition, for the and modes where 
the event samples are sizable, the measurements are repeated in bins of the 6-hadron 
kinematic variables. None of these checks reveals a systematic bias. 

5 Conclusions 

The 6-hadron masses are measured using data collected in 2010 at a centre-of-mass energy 
of y/s = 7 TeV. The results are 

M(fi+ ^ J/ipK+) = 5279.38 ± 0.11 (stat) ± 0.33 (syst) MeV/c^ , 

M(5° ^J/^pK^*'>^) = 5279.58 ±0.15 (stat) ±0.28 (syst) MeV/c^ , 

M(5° ^ J/^Ij(I)) = 5366.90 ± 0.28 (stat) ± 0.23 (syst) MeV/c^ , 

M(AO -^J/i)A) = 5619.19 ±0.70 (stat) ±0.30 (syst) MeV/c^ , 

where the B^ result is obtained as a weighted average of M{B^ — )■ J/tpK*^) = 5279.58 ± 
0.17 ± 0.27MeV/c2 and M{B° J li^Kl) = 5279.58 ± 0.29 ± 0.33MeV/c2 assuming all 
systematic uncertainties to be correlated, except those related to the mass model. The 
dominant systematic uncertainty is related to the knowledge of the average momentum 
scale of the tracking system. It largely cancels in the mass differences. We obtain 

M{B^ J/^fs:(*)0) -M{B+ J/i)K+) = 0.20 ± 0.17 (stat) ± 0.11 (syst) MeV/c^ , 
M(50 ^ J/^(f)) -m\b+ J/^K+) = 87.52 ± 0.30 (stat) ± 0.12 (syst) MeV/c^ , 
M(AO ^ J/^A) -M(5+ ^ J/^K+) = 339.81 ± 0.71 (stat) ± 0.09 (syst) MeV/c^ , 

where the 5° result is a combination of M(5° ^ J/^K*°) - M{B+ J/ipK^) = 0.20 ± 
0.20±0.12MeV/c2 and 1/(5° -> J/ipK^)- M{B+ J/^K+) = 0.20±0.31±0.10 MeV/c^ 
under the same hypothesis as above. 

As shown in Table 4, our measurements are in agreement with previous measure- 
ments [ , !]• Besides the difference between the B~^ and 5° masses they are the most 
accurate to date, with significantly improved precision over previous measurements in the 
case of the B^ and A° masses. 
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Table 2: Systematic uncertainties (in MeV/c^) on the mass measurements. 



Source of uncertainty 






5°-^ 








J/^K+ 


J/il^K*^ 


J/il^Kl 




Jli)k 


Mass fitting: 












- Background model 


0.04 


0.03 


<0.01 


0.01 


<0.01 


- Resolution model 


0.01 


0.02 


0.06 


0.02 


0.07 


- Radiative corrections 


0.01 


0.01 








Momentum calibration: 












- Average momentum scale 


0.30 


0.27 


0.30 


0.22 


0.27 


- Tj dependence of momentum scale 


0.04 


<0.01 


0.09 


0.03 


0.02 


Detector description: 












- Energy loss correction 


0.10 


<0.01 


0.05 


0.03 


0.09 


Detector alignment: 












- Vertex detector (track slopes) 


0.05 


0.04 


0.04 


0.03 


0.04 


Quadratic sum 


0.33 


0.27 


0.33 


0.23 


0.30 



Table 3: Systematic uncertainties (in MeV/c^) on the differences of mass measurements, ex- 
pressed with respect to the J/'^K^ mass {e.g. the last column gives the systematic 
uncertainties on M(Ag J/ipA) - M{B+ J/il;K+)). 



Source of uncertainty 


50 ^ 










J/iPK*^ 


J/^Kl 






Mass fitting: 










- Background model 


0.05 


0.04 


0.04 


0.04 


- Resolution model 


0.02 


0.06 


0.02 


0.07 


- Radiative corrections 


<0.01 


0.01 


0.01 


0.01 


Momentum calibration: 










- Average momentum scale 


0.03 


<0.01 


0.08 


0.03 


- r] dependence of momentum scale 


0.04 


0.05 


0.01 


0.02 


Detector description: 










- Energy loss correction 


0.10 


0.05 


0.07 


0.01 


Detector alignment: 










- Vertex detector (track slopes) 


0.01 


0.01 


0.02 


0.01 


Quadratic sum 


0.12 


0.10 


0.12 


0.09 
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Table 4: LHCb measurements, compared to both the best previous measurements and the 
results of a global fit to available &-hadron mass data [ ']. The quoted errors include statistical 
and systematic uncertainties. All values are in MeV/c^. 



Quantity 


LHCb 

measurement 


Best previous 
measurement 


PDG fit 


M{B+) 

M(AO) 
M(50) - M{B+) 
mIb{) - M{B+) 
M(AO) - M{B+) 


5279.38 ± 0.35 
5279.58 ± 0.32 
5366.90 ± 0.36 
5619.19 ±0.76 
0.20 ±0.20 
87.52 ±0.32 
339.81 ±0.72 


5279.10 ±0.55 [.] 
5279.63 ±0.62 [1] 
5366.01 ±0.80 [1] 
5619.7 ±1.7 [1] 
0.33 ± 0.06 [15] 


5279.17 ±0.29 
5279.50 ± 0.30 
5366.3 ±0.6 

0.33 ±0.06 
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